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We consider exact and quasi-exact solvability of the one-dimensional Fokker-Planck equation 
based on the connection between the Fokker-Planck equation and the Schrodinger equation. A 
, unified consideration of these two types of solvability is given from the viewpoint of prepotential 

■ together with Bethe ansatz equations. Quasi-exactly solvable Fokker-Planck equations related to the 

' s/(2)-based systems in Turbiner's classification are listed. We also present one s/(2)-based example 

^SJ , which is not listed in Turbiner's scheme. 
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' Without doubt the Fokker-Planck (FP) equation is one of the basic tools used to deal with fluctuations in various 
^ ' kinds of systems [l|. It has found applications in such diverse areas as physics, astrophysics, chemistry, biology, 
finance, etc. Owing to its wide applicability, various methods of finding exact and approximate solutions of the FP 
equations have been developed. 

One of the methods of solving the FP equation is to transform the FP equation into a Schrodingcr-likc equation, 
^ ■ and then solve the eigenvalue problem of the latter. This method, to be called the method of eigenfunction expansion, 
C/3 , is useful when the associated Schrodinger equation is exactly solvable; for example with infinite square well, harmonic 
oscillator potentials, etc. Several FP equations have been exactly solved in this way [ij. 

Unfortunately exactly solvable Schrodinger equations are rather limited in number, and hence the method of 
eigenfunction expansion is quite restricted. However, in real calculations more often than not a knowledge of a 
large but finite number of eigenfunctions is sufficient for a good approximate result. This being so, it would render 
^ the method of eigenfunction expansion more useful. Recently, in non-rclativistic quantum mechanics a new class of 
Q ' potentials which are intermediate to exactly solvable ones and non-solvable ones has been found. These are called 
O , quasi-exactly solvable (QES) problems for which it is possible to determine algebraically a part of the spectrum 
■ (eigenvalues and eigenfunctions) but not the whole spectrum 0, [3, IE @ • These QES Schrodinger equations thus 
^vq , extend the range of applicability of the eigenfunction expansion. 

^ ■ In this paper, we shall classify s/(2)-based QES FP equations according to their associated QES Schrodinger 
\ equations listed in and present a case which is not listed in We first discuss the connection between the 
• FP and Schrodinger equations in Section II. In Section III the basics of the exact and quasi-exact solvability of 
25 . Schrodinger equations are reviewed from the viewpoint of prepotential together with Bethe ansatz equations, and the 
way to construct the corresponding QES FP equations is presented. Some examples of QES FP equations constructed 
' in this manner are discussed in Section IV. Section V concludes the paper. 
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II. FOKKER-PLANCK AND SCHRODINGER EQUATIONS 

In one dimension, the FP equation of the probability density P{x, t) is jy 

^P(x,t) = £P(x,t), 

The functions D'''^\x) and D^'^\x) in the FP operator C are, respectively, the drift and the diffusion coefficient (we 
consider only time-independent case). The drift coefficient represents the external force acting on the particle, while 
the diffusion coefficient accounts for the effect of fiuctuation. The drift coefficient is usually expressed in terms of 
a drift potential $(a;) according to D'^^^x) = — ^'(a;), where the prime denotes the derivative with respect to x. 
Without loss of generality, in what follows we shall take D^"^^ = 1. The stationary solution of the FP equation is 
Po(x) = exp(-$(a;)), with / Pii{x) dx = 1. 
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The FP equation is closely related to the Schrodinger equation jT|. To see this, let us define tp{x,t) = e^/^P{x,t). 
Substituting this into the FP equation, we find that ij: satisfies the Schrodinger-like equation: 




V'. (2) 

As we shall see in what follows, the function $/2 plays a fundamental role in our discussions. So let it be denoted 
by W{x) = $(a;)/2. Following 3, we shall call W{x) the prepotential. Letting ip{x,t) = exp(— A<)(/)(a;) and (/)o = 
exp{—W), we obtain 



= A0, (3) 



where 



--S + n^), (4) 
V{x) = W'ixf -W"{x). 

Thus (/) satisfies the time-independent Schrodinger equation with Hamiltonian H and eigenvalue A, and 00 is the zero 
mode of H: H4>q = 0. Since only derivatives of W{x) appear in V{x), W{x) is defined only up to an additive constant. 
We choose the constant in such a way as to normalize 0o(2;) properly, J /p^ixY dx = 1. For simplicity of presentation, 
we consider the cases in which the ground state wave functions are square integrable, that is the corresponding FP 
operators have the normalizable stationary distribution. 

Now comes the essence of the method of eigenfunction expansion for solving FP equations. If all the eigenfunctions 
(j)n {n — 0, 1,2, . . .) of H with eigenvalues A„ are solved, then the eigenfunctions P„(x) of C corresponding to the 
eigenvalue — A„ is P„(a;) = 4)Q{x)(j)n(x). The stationary distribution is Pq = (jx^ = exp(— 2VF), which is obviously 
non-negative, and is the zero mode of C: CPq = 0. Any positive definite initial probability density P{x,0) can be 
expanded as P{x, 0) = 0o(a;) J2n Cn(t>nix), with constant coefficients c„ (n = 0, 1, . . .) 

<j)„{x) {(b^\x)Pix,0)) dx. (5) 

Then at any later time t, the solution of the FP equation is P{x, t) — (j)o{x) ^„ c„0„(x) exp(— A„t). 

It is now evident that the exact and quasi-exact solvability of FP equations are directly related to the exact and 
quasi-exact solvability of the corresponding Schrodinger equations. 

We shall classify the types of exactly solvable and QES FP equations associated with Schrodinger equations which 
are factorizable (or supersymmetric) . For this, it is desirable to have a unified approach to deal with both kinds of 
solvability. It appears to us that the approach based on the so-called prepotential serves this purpose well. We shall 
therefore first review exact and quasi-exact solvability from this viewpoint below. 



III. EXACT AND QUASI-EXACT SOLVABILITY FROM THE POINT OF VIEW OF PREPOTENTIAL 

WITH BETHE ANSATZ EQUATIONS 

Let us now briefly review the essence of quasi-exact and exact solvability of Schrodinger equations. We choose 
to proceed from the point of view of prepotentials together with the Bethe ansatz equations [7] (see also Q). The 
merit of this approach is that one needs not to first assume the form of the potential of the system concerned, as all 
information about the system is contained in the prepotential and the solutions, or roots, of the Bethe ansatz equations 
i. Also, in this approach exact and quasi-exact solvability can be treated on the same footing. Furthermore, such 
approach facilitates extension of the QES theory from the Schrodinger equation to equations for multi-component 
wave functions, such as the Pauli and the Dirac equation Q. 

Suppose (j)o = e~^°^^^ is the ground state, with zero energy, of a Hamiltonian Hq: i?o0o = 0. This is the simplest 
example of quasi-exact solvability. This implies that the potential Vq is completely determined by Wq: Vq — Wq^ — Wq , 
and consequently, the Hamiltonian is factorizable: 
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This fact can also be considered as the very base of supersymmetric quantum mechanics. Consider now a wave 
function (N: positive integer) which is related to (f>o of Hq by (f>pf = 4>o(f>N- 

(pN (z - zi){z - Z2) ■ ■ ■ {z - zn). (7) 

Here z = z{x) is some function of x. The function cj)]^ is a polynomial in an (A'^ + l)-dimensional Hilbert space with 
the basis (1, z, z^, . . . , z^). One can rewrite (pN as 

(/>Ar = exp {-Wn{x, {zk})) , (8) 

with 

N 

Wn{x, {zk}} = Wo{x) - ^ In \z{x) - Zk\. (9) 

k=l 

Now we form Wj^ — W^. This results in a Riccati equation 

W'j^ix, {zk}) - W'j^ix, {zk}) = VNix, {zk}) - XN{{zk}), (10) 

where Vjv is in general a function of N, x and {zk}, and Ajv is a real constant depending on N and {zk} (dependence 
of V^v and Aat on other parameters of Hq is understood). It should be noted that the definition of Vn and Ajv has 
ambiguity, namely, they can both be shifted by an arbitrary real constant a: 

Vn ^ Vn + a, XN^^N+ct- (11) 



This corresponds to the redefinition of the zero point of Vn- For the theory we are considering, the r.h.s. of ([Tt 
is generally a meromorphic function of z with at most simple poles. Let us demand that the residues of the simple 
poles, Zfc, fc = 1, . . . , iV should all vanish. This will result in a set of algebraic equations which the parameters {zk} 
must satisfy. These equations are called the Bethe ansatz equations for {zk}- For each value of N {N = 0, 1, 2, . . .), 
there are N + 1 possible sets of solution {zk} of the Bethe ansatz equations 

With {zk} satisfying the Bethe ansatz equations, the r.h.s. of pU|) will have no simple poles at {zk} but it 
still generally depends on {zk}- Suppose now that the N + 1 sets of roots give the same form of the potential 
Vn{x, {zk}) — Vn{x) in the r.h.s. of ([10]), but with + 1 values of Ajv. Then the Riccati equation (fTO|) implies that 
we have a new eigenvalue problem with Hamiltonian TIn' 

'Hn4>n = Xn4>n- (12) 

The 1 sets of {zk} define iV+ 1 eigenfunctions (f>N of TYat with eigenvalue Ajv. In other words, iV+ 1 eigenfunctions 
and eigenvalues of the Hamiltonian Tijv are solvable algebraically, and Tijv is said to be QES. If the wave function 
(/>j^'' = exp(— W^^) corresponding to the lowest value Xn = Aj^"* has no node, then it is the ground state of T^at, and 
is annihilated by the Hamiltonian Hm- 



H^^Hj,- Ai,") = + - <^", (13) 



Hn<I>^n^ = 0. (14) 

Now if Vn = V is not only independent of the solutions of the Bethe ansatz equations {zk}, but is also independent 
of N, then the Hamiltonian H is exactly solvable, as solvable states can be found for any integral value of N. 

From the Schrodinger equation TLniPn = X^fpN, it is also seen that the polynomial part (jjN satisfies the equation 
hN(t>N = Ajv^AT, where = 4>q'^'Hn4>o^ and 0o is the ground state of Hq corresponding to TV = 0. From the previous 
discussions, we know that /ijv has an algebraic sector with 7V+ 1 eigenvalues and eigenfunctions, which are polynomials 
in an (iV + l)-dimensional Hilbert space with the basis (1, z, z^, . . . , z^). The quasi-exact solvability of Hm is said to 
have a Lie-algebraic origin, if can be expressed as 

hN = Y^ CabJ^J^ + ^^■^'^ + constant, (15) 



where Cab, Ca are constant coefficients, and the J° are the generators of some Lie algebra. QES Hamiltonians based 
on the Lie-algebra sl{2) have been classified in [1, 0]- 
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There exists a different situation when the N + 1 sets of roots {zk} are plugged into the l.h.s. of ITOl) . Instead of 
resulting in the r.h.s. of pUj) a single potential Vn with + 1 values of Ajv, N + 1 potentials with a single value 
of Xn are obtained, when the zero point of the eigenvalue is properly adjusted. These are QES systems for which a 
set of iV + 1 potentials differing by the values of parameters have the same eigenvalue of the j-th eigenstate in the 
j-th potential (j = 1, 2, . . . , + 1). Such kind of QES systems were termed the second type QES systems, and those 
discussed previously were called the first type in 3]. For the second type QES quantum systems which have sl(2) 
Lie-algebraic origin, it is (h^ — E) / p{x) instead of ftjv, where p{x) is some x-dependent scale factor, that is expressible 
in terms of certain quadratic combination of the generators J^'s. In this sense, one can say that for the first type 
systems, the scale factors p are simply some real constants. 

Combining the discussions in this section and those in the preceding section, we arrive at a recipe of constructing 
QES FP equation from a known QES Schrodinger equation. We shall be interested only in the first type QES systems, 
since for the second type QES systems the corresponding Wn only define N + 1 FP systems with + 1 different 
drift potentials = 2Wn having the same eigenvalue. While these are QES FP systems, it is not useful for good 
approximate calculations using the method of eigenfunction expansion as only one eigenstate is known. 

Suppose a prepotential Wq and a corresponding set of Bethe ansatz equations define a first type QES Schrodinger 
equation in the sense discussed above. Then the drift potential $o = 2Wo defines a QES FP equation with only 
one exactly solved state 0o, and the stationary solution of the FP equation is Po{x) — (j>Q. To obtain QES FP 
equations admitting a larger number of solvable states, wc construct a new prepotential Wn from Wq and the N 
roots {zi, Z2, ■ ■ ■ , Z]\[} of the Bethe ansatz equations according to The set of roots {z^} to be chosen is the set for 
which the value of eigenvalue \m is the lowest. If the corresponding wave function 0^-* has no node, one can use it to 
define a corresponding FP operator: 

Cn^-^^n'^Hn^n^-'. (16) 

Then Cn defines a QES FP equation having 1 solvable states, with a drift potential = 2W^''. As mentioned 
before, if is independent of N, then is exactly solvable, and so is the corresponding FP equation. 
We shall illustrate this procedure by some examples below. 

IV. EXAMPLES 
A. Exactly solvable cases 

It is clear that FP equations transformable to exactly solvable Schrodinger equations can be exactly solved. Par- 
ticularly, all the shape- invariant potentials in supersymmetric quantum mechanics, as listed in Table 4.1 of [lo| . or 
equivalently. Table 5.1 in [ll[, give the corresponding exactly solvable FP systems. One needs only to link the pre- 
potential W [m in the Schrodinger system with the drift potential or drift cocfhcient D'-^\ in the FP system 
according to the relation $ — 2W and D'-^-' = — We shall not list all of them here, but simply mention that the 
FP equation for the Ornstein-Uhlenbeck process corresponds to the shifted oscillator potential in quantum mechanics, 
and that for the Rayleigh process corresponds to the three-dimensional oscillator potential. 

B. sZ(2)-based QES cases listed in Turbiner's classification 

Similarly, one can identify all QES FP equations which are based on the sl{2) Lie algebra from the list of ten classes 
depicted in Table 1 of Q. Of the ten classes, classes II, III, V, VIII and IX belong to the second type QES systems 
which, as discussed previously, are not useful as far as the method of eigenfunction expansion is concerned. Class 
X gives a periodic potential, which is also not of concern in this paper, as we are interested in normalizable wave 
functions so that the distribution density of the FP equation is also normalizable. 

So we are left with only classes I, IV, VI and VII in the list of Q which give QES FP systems having iV -I- 1 
eigenstates with arbitrary N. To obtain these four QES FP systems, one needs only to identify Wq with the gauge 
function g of each class listed in [s*]. This gives a QES FP system with only one solvable eigenstate. For QES FP 
systems in the same class but with higher values of N, wc construct Wn from Wq and the roots of the Bethe ansatz 
equations {zfc} (with the function z = z{x) also from [3]) according to the procedure in 

Let us illustrate this by class VII QES FP system. It turns out that this case is just the one-particle case of the 
so-called rational BC type Inozemtsev model Q . 
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In this case z{x) = x"^, and the prepotential ([9]) is given by 

Wn{x) ^Wo + W, (17) 
Wo = ^a;^ + ^x^ - jinx + Cn, x > 0, a > 0, 7 > 0, (18) 

N 

W = -Y,^n\x'' - Zk\. (19) 

fe=i 

Here Cn is a constant term necessary for the normahzation of the ground state wave function (l)^^\x). By plugging in 
all + 1 possible solutions {zk} into Vn in ([10]), we find that Vn is independent of {^fc}. The Bethe ansatz equations 
removing the simple poles at {zk} read 



Zk 

7.1. 

l^k 

and Ajv in terms of the roots Zk 's is 

N 

Zk 



2oz^ + 26zfc-(27+l)-4^ ^ 0, fc = l,...,7V, (20) 



^ 1 

Ajv = 2(27 + l)^-, (21) 



fc=i 

or equivalently, using the Bethe ansatz equation, 

N 



Aat = 46/^ + 4a^Zfe. (22) 

k=l 

Among all the solutions, we choose those {zf^} which are all negative 

Zk<0, k^l,...,N. (23) 
This corresponds to the ground state, since the eigenvalue A^^' is the lowest, and the wave function 

N 

(0) ^^-Cix^+^x')^! TTc^2 



"AT 



l[{x'-Zk) (24) 



k=l 



has no node. The fact that the choice ((23)l is possible can be understood as follows. Recall that the number of QES 
states is + 1 for each N. Since the factor nfc=i(^^ ~ ^fc) wave function is even in x, there can be at 

most N zeros in the half line a: > 0. According to a well known theorem in quantum mechanics, no two eigenstates 
can have the same number of nodes. This implies that the number of nodes of the + 1 QES states here must be 
0, 1, 2, . . . , TV. Hence the ground state is in the QES sector, and the choice (j23p is guaranteed. A similar argument 
can be applied to the other cases. 

The corresponding QES FP equation is then defined by the drift potential 

N 

$7v = -X* + bx^ - 27lnx - 2^ In \x'^ - Zk\+ constant. (25) 

k=l 

The constant is chosen so that the stationary distribution of this FP equation 0^''^ = exp(— $Ar) is normalized. We 
have thus obtained a QES FP equation corresponding to the class VH s^(2)-based QES Schrodinger equation in 

QES FP equations of the other three classes can be treated accordingly. The main defining information is given in 
Table 1. 



C. s?(2)-based QES case not listed in Turbiner's classification 

We now present a QES case which is also based on the sl{2) algebra, but not listed in [sj. This is related to the 
trigonometric BC type Inozemtsev model with one degree of freedom |7| . 
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TABLE I: Three other QES FP systems associated with the sZ(2)-based QES Schrodinger equations in 3]. The function z{x), 
prepotential Wo{x) (—00 < a; < 00), Bethe ansatz equations, and the eigenvalues Ajv are listed. The drift potential defining 
the QES FP equation with A'' + 1 solvable states is <l>jv = 2Wn, where the prepotential Wn is constructed from ^ with Wo 
and the roots of the Bethe ansatz equations. All parameters are real. 

Class z{x) Wo Bethe ansatz equations (fc = 1, 2 . . . , TV) Ajv 



azl 



(a > 0, a,c > 0, V 6) 



Zk ~ C 

= 



IV 



cosh 



■ cosh 2ctx + — In cosh ax 

a 

(a > 0, c > 0, V a) 



(a + |a)zi + (c- 



-2a E, 



Ol)zk ■ 





2ac- 



VI 



4 ^ 2 

(a > 0, V & or a > 0, b > 0) 



az'i + bzk - ^ 
-2 V ^fc 



= 



The prepotential in this case reads 



Wn — ^ cos 2x — ^ In I cot x\ ~ c log | sin x 

N 



-^Inlsin^x-Zfcl +CAr, < x < -. (26) 
fc=i 

All the parameters a, b and c are real and satisfy a > 0, c > | > 0. Again, Cn is a constant term necessary for the 

normalization of the ground state wave function (j>^^\x). The first three terms in (j26p are just Wq which gives only 
the ground state solution. The parameters {zk} are determined by the following Bethe ansatz equations: 

(4azfe + 2c)(l - Zk)-b+l-2zk + 4zfe(l - Zfe) V = 0, 

fc = l,...,iV. (27) 
With these solutions, the potential Vn and eigenvalue Aat in the Riccati equation PU)) are 

Vn{., {z,}} = - ^ - SaiVsin^ x, (28) 

and 

N 

XN{{zk})=SaY,Zk+iN{N + c). (29) 
fc=i 

Again we choose those {2^:} which are all negative 

Zfe<0, k = l,...,N. 
This corresponds to the ground state and its wave function has no node: 



N 



(X e~'^''°''^''{cotx)^smxy Y[{sm^ x ~ Zk). (30) 



k=l 



The corresponding QES FP equation is then defined by the drift potential $jv = SWj^'. 
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The prepotential Wn in (|26p . and hence the potential Vn in ((28)) are not hsted in the sU2)-based classification 
scheme in [j]. Nevertheless, it turns out that this system is also related to sl{2) Lie algebra This can be seen as 
follows. Using (/)o = exp(— Wq) and Vat, we obtain the transformed Hamiltonian Hn = 4'o^T^n4'o as 



I AT 



dWo d „ . 2 
dx-^ dx dx 

d^ f b \ d , 

---^ - 2 asin2a; — — + ccota: SaTVsin^ x. (31) 

dx'' \ sm2a; / dx 



In terms of the new variable z{x) — sin^ x, hfq becomes 



Hn = -4z(l - z)—r + 2 \Aaz^ + 2 (c - 2a + 1) z + 5 - 2c - ll — - SaNz. (32) 
dz-' dz 



Now it is easy to check that eq. ([32]) can in fact be rewritten in the form (fT5|) as 

= - AJ^J- + 8aJ+ + 4 (c - 2a + 1 + Af) J" + 2 (6 - 2c - 1 - iV) J" + 2A^ (c - 2a + 1 + TV) , (33) 

where the generators of the Lie algebra sl{2) are given by 

J+ = z^^-Nz , 
dz 

J- = A , (34) 

dz 

This proves that the trigonometric BC type Inozemtsev model is indeed a QES system based on sl{2). 



V. SUMMARY 



In summary, we have discussed exact and quasi-exact solvability of the FP equation based on the corresponding 
solvability of its associated Schrodinger equation. We give a unified treatment of these two types of solvability from 
the viewpoint of prepotential together with Bethe ansatz equations. Examples of QES FP equations related to the 
s/(2)-based systems listed in jSj are classified, and one example which is also based on the sl{2) algebra but not listed in 
3 is also presented. We note here that, as shown in 0, the concepts of quasi-exact solvability and the so-called higher 
derivative, or nonlinear, or M-fold supersymmetry |14| are equivalent. Hence, all the QES FP equations discussed 
here also possess such kind of extended supersymmetry. 

As mentioned before, we only discuss those QES FP systems which can admit a large number of QES states so 
that they are useful for consideration when approximate calculations are required. Hence, we have not considered 
the so-called second type QES systems defined in as they admit only a single solvable state and are therefore of 
less interest in practical calculations. Nevertheless, these QES quantum systems still give the corresponding QES FP 
equations. We hope that in the future more non-s^(2) QES FP systems can be found. 
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